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A  TURBULENT'  MHD  BOUNDARY  LAYER 
IN  A  TRANSVERSE  MAGNETIC  FIELD 

I.  P.  Ginzburg  and  L.  I.  Skurin 

In  this  work  we  examine  the  turbulent  boundary  layer  of  a 
conducting  fluid  on  a  dielectric  plate  in  a  transverse  magnetic 
field.  We  investigate  both  nongradient  and  gradient  flow.  The 
basis  of  the  solutions  is  the  semiempirical  turbulence  theory,  and 
we  use  the  calculating  method  shown  in  [1]  for  a  ordinary  boundary 
layer. 

Considering  that  the  magnetic  Reynolds  numbers  are  small, 
we  can  write  the  equations  of  the  boundary  layer  in  crossed  electrical 
and  magnetic  fields  associated  with  the  plate  in  the  form 


■fc  +  T?-.*  ■ 

(1) 

(2) 

(3) 

fWj,  gjr  4*  PVy  *=  *1'  ^vx)  4*  (E*  vx&)> 

(4) 

•  »  p*=  fRT, 

(5) 

where  E„  is  the  external  electrical  field  strength;  B  is  magnetic 
induction;  c  is  conductivity  of  the  fluid;  q  is  the  transverse 
component  of  the  heat  flow  vector;  t  is  the  xy  component  of  the 
stress  tensor;  H  is  the  total  enthalpy  per  unit  of  mass;  the  remaining 
designations  are  those  commonly  used. 
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Considering  (3)  and  bearing  in  mind  that  outside  the  boundary 
3V 

layer  -=r~  ~  0  and  -r-  =  0,  from  (2)  we  get 
o  y  oy 

—  ~q*®£s*  (6) 

where  u  is  the  velocity  on  the  outer  boundary  of  the  boundary  layer, 

and  the  subscript  M0”  pertains  to  the  outer  ^boundary.  From  (4)  , 

which  pertains  to  external  flow,  considering  that  outside  the 

boundary  layer  q  *  x  ■  ~  «  0,  we  get 

®  y 

whence  it  follows  that  when  Ez  *  0  the  total  heat  content  outside 
the  boundary  layer 

Ht  —  c,T%  +7  “  const.  (7) 

Examination  of  (6)  leads  to  various  statements  of  the  problem. 

1.  Nongradient  flow  ~  3^  In  this  case  the  velocity  on 

the  outside  cboundary  is  determined  from  the  relationship 

M  -f  9|fltu  _  aJBE' = 0,  (8) 

where  the  external  magnetic  and  electrical  fields  should  be  given. 

In  the  particular  case  Oq  «  0,  from  (8)  it  follows  that  u  ®  const. 

2.  Gradient  flow  —  ^  ^  0.  In  this  case  the  pressure  gradient  , 
is  determined  from  (6)  with  given  o0,  u,  B,  and  Ez.  Let  us  examine 
certain  of  the  simplest  cases  of  the  indicated  statements  of  the 

,  problem.  When  solving,  the  values  of  t  and  q  which  enter  into  (2) 
and  (4)  will  be  determined,  starting  from  the  semiempirical  turbulence 
theory,  using  the  two-layer  scheme  of  the  boundary  layer.  According 
to  this  theory, 

1 + w,  -  r  [-pr  •  %  +  (1  -  w)  »*■ isf\ : 

where  6^  is  the  thickness  of  the  laminar  sublayer;  y,  X,  cp  are. 
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respectively,  the  coefficients  of  dynamic  viscosity,  thermal  conduc¬ 
tivity,  and  heat  capacity  at  constant  pressure;  T  is  temperature; 
cdu  k 

Pr  *  Pr«-  s  IT*  kt  are  universal  constants  of  the  semi- 

empiricai  turbulence  theory.  Henceforth  we  will  consider  Pr  and  cp 
constant  and  k  *  . 


It  is  assumed  t-fc*>t  when  crossing  the  boundary  of  the  laminar 

sublayer  the  values  x,  q,  h  *  e  T,  and  v  do  not  undergo  a  discon- 

P  X 

tinuity,  while 


r° 


(9.) 


where  k.^  is  a  universal  constant. 


For  an  approximate  solution  to  system  (l)-(5.)  we  set 

H—Bx  +  BtV' r-f  Btvl  when  y <  $i» 

//'■*  Cjt/^  when  y  ^  V.10) 


Coefficients  A^,  A2 ,  B^,  B2,  B^,  C^,  and  C2  are  determined  using 
the  boundary  conditions,  the  motion  and  energy  equations,  relationship 
(9),  and  the  condition  x  «  0  when  y  *  6  (6  is  the  boundary  layer 
thickness) . 

The  solution  for  such  a  scheme  for  a  nonconducting  fluid  is 
given  in  [1].  Henceforth  we  will  use  the  results  of  this  work. 


Nongradient  flow,  Ez  *  0,  0  =  const.  We  have  the  same  boundary 
conditions  as  for  the  boundary  layer  of  an  nonconducting  fluid: 

H-HW-C,TW ; 

ymml,  Vj,**  U%  H**  ' 


where,  according  to  (7),  HQ  *  const.  Let  us  define  u.  The  density 
on  the  outer  boundary,  according  (3),  (5),  and  (7),  can  be 


represented  in  the  form 


CpTw 

~w 


(12) 
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Substituting  this  expression  into  (8),  where  Ez  =*  0,  we  get  an  equation 

* 

for  determining  the  velocity  at  the  outer  boundary 


Jyr- 


~W' 


l~w; 


Substituting  this  expression  into  (u),  where  Twpw  =  const  (since 
p  *  const),  from  the  last  equation  we  get 

Arthu  —Arthu, 


P./57 J%Hm' 


(13) 


where  «= 


yifh'  VWt; 


3X 

form 


The  integral  pulse  relationship,  corresponding  to  (2),  where 

*  0,  E_  *  0,  with  consideration  of  (8),  can  be  written  in  the 
z 


where 


(14) 


(15) 


Since  in  this  case 


il 

3y 


0  when  y  *  0,  the  dependence  of  t  on  y/6 

can,  approximately,  be  represented  in  the  form  t  *  tw|l  -  ^ j  ,  arid 

then  for  the  velocity  profile  we  get  the  ;same  expression  as  in  (l); 
6**  on  c  =  yi  -H—  will  be  the  same  as  in  the  absence  of  a  magnetic 

"Vpw 

field. 


Calculation  of  the  integral  on  the  right  side  of  (14),  using 
the  velocity  profile  obtained  in  [1],  gives 


j(  i  -*)«--*?  + 


+'■ 


~^-r—  (rfjcosdt+^n^*)  . 
■j*r  ‘  . 


■where  d,  *  arcsin  — 
1  a 


1  - 


2u  ' 


,  d. 


arcsin  -~o>  and>  according  to  [1], 
2au^ 


,  Dae£)(u,  Hw,  Pr),  a*=a(u,  7/*,  «*); 


/an;  _  _ 

aEa(;,  u,  ftw,  Pr),  a==a(C,  u,  //»,  Pr). 
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The  values  of  pA  and  p  as  functions  of  x  are  determined  by 
0  w 

relationships  (12),  (13),  (5),  and  (11).  Thus,  the  problem  reduces 
to  numerical  integration  of  (1*0. 


Let  us  examine  an  incompressible  fluid.  Pulse  relationship  (14) 
is  written  in  the  form 

(16) 

the  velocity  on  the  outer  boundary,  according  to  (13),  will  be 
defined  by  the  expression 

the  width  of  the  lost  pulse,  according  to  [1], 

(17) 

2 

Setting  in  this  last  relationship,  approximately,.  1  —  —•  • 

■  1  -  =0.8  and  substituting  (17)  into  (16),  we  get  the 

?Jcp 

following  equation  for  determining  5  as  a  function  of  x: 

2F  (“**)“  (l8) 

To  obtain  an  approximate  analytical  solution  to  (18)  let  us  replace, 
considering  Sx  small,  u  by  uQ  on  the  left  side  of  (18)  (the  error  in 
determining  ?  when  S  £  0.4  does  not  exceed  6%).  We  then  get  the 
solution  in  the  form 


where 


Re^l  -5,  +  y  , 

»c 

e  j  ijWij  =  e‘ « (K - 1  )*  +  ^  -f  2. 


(19) 


(20) 


Considering  kc  to  be  a  sufficiently  large  value,  let  us 
disregard,  in  this  last  equality,  the  second  and  third  terms  on  the 
right;  then,  using  the  approximation 

1).  «i-  const,  /i,  =  const,  (gi) 

we  get  .  ' 

{M) 

whence 
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:a?V,nJ 


A4  exp  {*i  -f  itt  —  n,  —  l) ' 

0,8*, 

_ J_c  _r  _  '  c 

%*-  ■ r  *,*  0jt* 


Re 


<•]- 


(23) 


For  values  of  kc  =  6-12,  n1  =  2*1$  and  n2  =  1.2-5.  Setting  k-j.  -  4.29 
and  k  *  0.39,  just  as  for  the  usual  turbulent  boundary  layer,  from 
the  last  relationship  we  get 

®  (j.o  —  2, OS,. 

i 

In  [2],  where  there  is  obtained  a  solution  for  the  examined  problem 
from  concepts  of  the  theory  of  dimensionality,  there  is  given  an 
analogous  relationship  with  the  coefficient  I.87. 


Haying  set  v* 


whence,  using  (23) 


we  have 


V®  v* 


we  find 


(1  -sx), 


0.434(1-5,) 


.  f^expf*!  -f /•»- 

'*[ —  w 


~  ReJf 


(2 1) 


Figure  1  gives  the  limiting  curve  constructed  for  the  case  S  ■  0,  on 
the  right  side  of  (24).  The  same  figure.'  shows,  for  comparison,  the 
corresponding  curve  obtained  in  [2], 


Fig.  1.  Rate  of  decrease  of  dy¬ 
namic  velocity  per  unit  of  parameter 

s4§£  ’ o,>  0  * const.  • 

1  -  from  data  in  [2]j  2  -  present 
work. 


Let  us  introduce  a.  formula  for  the  local  friction  coefficient: 

(25) 

r°o 

From  (21),  with  sufficiently  large  n,  we  have 
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*  1 


(26) 


which,  after  substitution  into  (25)  and  use  of  (22),  gives 

•2*t  f**wp{*i4-/4  — «t -l)Kw 

C/-\W'*,  I  5^1  J  X 

x  (l  _  sx  +isi)5'*,(l  -  S,)»  Rc?“‘ , 


whence 


In  [3],  the  value 


'fB=0 


=  1  —  2.69x  +...  was  found  for  the 


laminar  boundary  layer.  Comparison  of  the  ratios  of  local  friction 
coefficients  for  laminar  and  turbulent  boundary  layers,  given  in 
Fig.  2,  shows  that  the  magnetic  field  has  a  stronger  Influence  on 
the  laminar  layer. 


Fig.  2.  Comparison  of  the  ratios 
of  local  friction  coefficients 
for  laminar  (1)  and  turbulent  (2) 


boundary  layers 


= 

dx 


0 ,  a 


const 


Nongradient  flow  E„  -  0,  o  =  0.  As  in  [4],  we  use  for  the 

z 

conductivity,  the  following  integral  relationship: 


(27) 
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whence  o0  =  0,  and  then,  as  a  result  of  (8),  the  velocity  on  the  outer 
boundary  is  constant  (u  =  const).  The  boundary  conditions  are  given 
by  relationships  (11). 


The  integral  pulse  relationship  corresponding  to  equation  of 
motion  (2),  where  =  0,  E  =  0..  and  a  is  determined’  by  (27),  will 

d  X  2* 

have  the  form 


<«**  _  2*»  a,  JX- 

~4x  PvU  Po“*  ■ 


(28) 


The  velocity  profile  in  the  boundary  layer  in  this  case  as  well 
can  be  taken  approximately  as  in  [1]. 


The  width  of  the  lost  pulse,  according  to  x[l],  can  be 
represented  approximately  by  the  expression 


&•*  _  .  bz  .  ^s.  i  /  Hw  c 

•  N  -  V  l-«‘ 


Since  u  =  u 


Q  -  const,  then  TQ  is  also  constant  by  virtue  of  (7). 


Henceforth  we  will  consider  that  Hw  *  const  (or  Tw  *  const);  then 

p  =  const,  H  =  const,  and  u  *  const.  Having  substituted  this  last 
w  w  w 

formula  into  (28),  and  making  simple  transformations,  we  get 


g*K  ES 


l 


5TW  ’ 


(29) 


where  Dx> 


ID 

«**>( 


It 

—~r 


In  the  case  of  a  uniform  magnetic  field  (B  =  const).,  equation 
(29)  reduces  to  the  quadrature 


where  S 


k2 

o  b  v 
W 

2* 

pwu0 


f 


ND^t'  +  l 


51'., 

Dx 


To  obtain  a  solution  in  elementary  functions,  we  can  use 
approximation  (21). 
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Let  us  examine  in  more  detail  for  the  possibility  of  comparison 
with  the  results  obtained  for  the-  laminar  layer,  the  case  when 
B  =  Bq//x,  Bq  *  const.  Let  us  set,  in  (29), 


where  S 


gwB0 

w 


T 

BP~  'Di(i-S«*'-1) 


!  * 
*  (aktyffK 


We  get 


or 


d  c**  1 


Re, 


whence , 


by  virtue  of  (20)  and  (21),  we  get 


(«*)**•*- 


This  last  relationship  makes  it  possible  to  transform  the  expression 
2  BQ 

for  B  to  the  form  B  =  and  is  the  solution  of  (29)  for  the  given 
distribution  of  the  magnetic  field. 


Using  approximation  (21)  we  get 


Di  (i-  &*•”*); 


X  n 

From  this  it  follows  that  when  S  ^  ex~  2,  there  is  separation  of 
the  boundary  layer. 


(30) 


Relationship  (26),  using  (30),  gives  a.  formula  for  the  local 
friction  coefficient: 


Cf  *5i  2i»  _  2  xa'Al _ fef.  .  V 

7  fu*  CT  Ht  [i-St”'-*) 


l-l 

*« 


2  2 

where  A  =  2k  a  ngexp 


f«o 

n2  +  r\1  +  1^  1-  1 


(3D 


n. 


D, 


n2 


From  this  the  coefficient 


of  friction  of  a  plate  of  length  L 


L  j _ ,  ,_j_ 

c/=  j* cfdx  -  A'RqI‘~  (1  —  Se”>~') 


(32) 
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Using  (27)  and  (29)  it  is  also  easy  to  calculate  the  coefficient 
bf  the  stresses  caused  by  three-dimensional  magnetic  forces  in  the 
ooundary  layer: 

i  •  1 _ | 

CM  — tt  f  f  oB'vJydx  —  *Re£  — — - T  • 

'  (i 

The  coefficient  of  resistance  of  a  plate  of  length  L 

.  _i_  ,  j_ 

Cc~C/+CM  m,  A Re^  ~  (1  -Stf**-')"*'  • 

The  heat  flow  on  the  plate  is  given,  according  to  [1],  by  the 
expression 


1  —  Pr 


(33) 


where  u.  =  ^  «,  ^ •  Substituting  (31)  into  (33)  ,  we  get  the 

distribution  of  the  heat  flow  along  the  plate. 


In  the  case  Pr  =  1,  from  (33)  it  is  easy  to  determine  the  heat 
flow  of  a  plate  of  length  L 

L 


Qw  *  X  f  9lPlPf> 


1 


whence,  considering  (32), 


r  ,_L 


B-0 


Figure  3  shows  the  dependence  of  the  coefficients  C^,  C^,  and 

C  on  parameter  S.  Figure  4  compares  the  heat  flow  ratios  with  and 
c 

without  a  magnetic  field  for  laminar  [4]  and  turbulent  boundary 
layers.  As  follows  from  Fig.  4,  the  magnetic  field  exerts  less 
influence  on  the  turbulent  layer  than  on  the  laminar  layer. 
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Fig.  *J.  Comparison  of  heat  flow 
for  laminar  (1)  and  turbulent 


(2)  boundary  layers  ^  *  0, 

Bq.1 

°0  =  B 


Gradient  flow,  a  *  const.  In  this  case  the  pressure  gradient 

is  determined  from  (6)  by  giving  the  values  u,  B,  and  E  .  By  virtue 

z 

of  the  constancy  of  the  pressure  across  the  boundary  layer,  we  can, 
substitute  (6)  into  (2);  then  we  get 


When  there  is  an  external  electrical  field,  it  is  more  difficult 
to  determine  the  parameters  of  the  compressible  fluid  on  the  outer 
boundary  since  the  condition  of  constant  total  heat  content  is  not 
satisfied.  If  we  consider  that  there  is  no  external  electrical  field 
the  temperature  and  density  on  the  outer  boundary  are  determined  from 
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the  condition  of  constant  heat  content  (7)  and  the  Clapeyron  equation 
(5).  From  (7)  we  have 


r  1 


<jr 


«  do 
c,  rfx  ’ 


(35) 


differentiating  (5)’  with  respect  to  x  and  using  the  last  relationships, 
we  get 

zr  *  r,  n  ( — t)  “  <7  s 

or,  replacing  ^  from  (6)  (where  Ez  =  0),. 

Ut- 06) 

Thus.,-  with  given  u  (and  also  aQ  and  B),  the  temperature  and  density 
on  the  outer  boundary  are  determined,  respectively,  by  the  first 
of  relationships  (35)  and  by  relationship  (36). 


In  the  case  of  an  incompressible  fluid,  if  the  conductivity 
does  not  depend  on  the  temperature,  the  dynamic  and  thermal  problems 
are  separate.  In  particular ,.  when  a  «  const,  the  electrical  field 
does  not,  according  to  (34),  influence  flow  in  the  boundary  layer  and 
enters  only  into  (6)  which  connects  the  pressure  .gradient  and  the 
velocity  on  the  outer  boundary.  Let  us  examine  this  case  in  detail. 

Coefficients  and  k^  of  approximation  polynomial  (10)  for 
friction  are  determined  from  the  conditions 

y-8.,-0,g-0 

(the  last  equality  follows  from  (34));  then 

*  =  (37) 

Starting  with  this  last  relationship  and  using  the  propositions  of 
the  semiempirical  turbulence  theory  stated  previously,  we  find 

-KV 

y>8ii - f); 

«*  1  *! . 
v  Tv* 

Joining  the  velocity  profiles  on  the  boundary  of  the  laminar 
sublayer  we  get 
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(38) 


The  integral  pulse  relationship*  corresponding  to  (34),  is 
written  in  the  form 


a*  *  v  t-  a,  g(t** -}•  I*)  til* 

a  "TReJ" Re,  ‘  *  V«*  *  * 

U0X 

where  Re  =  — 5**  and  5*  are  determined  from  (15)  (pq  *  p)Ji 


(39) 


In  view  of  the  small  thickness  of  the  laminar  sublayer,  let  us 
calculate  approximately  the  integrals  in  (15),  using:  only  the 
velocity  profile  corresponding  to  the  turbulent  part  of  the  boundary 
layer.  Having  integrated,  and  considering  (38),  we  get 

£  -  —  -  &*'-**“  ('  -  «)•  («0) 

Let  us  calculate  the  constant  coefficients  in  (38)  and  (39),  and  the 
factor  in  parentheses  in  (40),  which  is  the  ratio  d**/i*»  let  us 
average,  in  the  proposed  range  of  change,  c  and  designate  it  by  ic 
(in  the  first  approximation  we  set  k  *  1).  We  will  have 

=  (41) 


Substituting  this  into  (39),  we  get,  finally,  the  equation 


,K»  **  _  0.29  ut  du  /  *  |  „x 


* 


(42) 


integration  of  which  gives  the  solution  of  the  problem  for  given 
u ( x ) ,  B(x ) . 


Let  us  examine  certain  analytical  solutions  of  (42).  When 
u  =  uQ/l  +  b  Rex  and  B  «  BQ/1  +  b  Rex  (b  >_  0,  uQ  *  const,  BQ  *  const), 
equation  (42)  reduces  to  the  quadrature 

K 


J  l+3,46^(*+I)-,^]nV>  ir,n<l+dRe-v)' 

°Bqv 

where  N  =  - 2* 

puQ 


(43) 


When  N  >  b(i<  +  1),  integral  (43)  diverges  at  the  point  defined 
by  the  relationship 

<*>•«“— srrJf+TF-  («> 
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■f 


Prom  this  it  follows  that  with  sufficient  large  N  the  parameters  of 
the  boundary  layer  will  have  limiting  values  that  are  defined  by  (44) 
When  N  £  b(.<  +  1).,  the  nature  of  the  change  in  parameters  will  be 
the  same  as  in  an  ordinary  boundary  layer. 


Let  us  examine  in  more  detail  the  case  b  «  0  Ou  *  Uq  *  const, ; 
B  *  Bq  »  const).  Let  us  set  k  *  0.8  (the  proposed  range  of  change 
of  kc  is  from  6  to  12).  Relationship  (43)  is  rewritten  in  the  form 


Ff 


fe--0'861^ 


Integrating  using  approximation  (21)  and  making  the  obvious  transfor¬ 
mations,  we  get 

jMC—  ft  Rt-V 

Nf11-'  '•  (45) 


whence,  as  a  result  of  (26), 


(46) 


W 

The  limiting  values  of  — *■  as  functions  of  N,  calculated  from  (44) 


(when  b  *  0),  are  given  in  Fig.  5.  (Naturally,  these  same  values 

can  be  obtained  approximately  from  (46)  when  Rex  ■+■  ®.)  This  same 

figure  gives  the  limiting  values  of  the  dimensionless  thickness  of 

Up.  6  —  v  fi  hO 

the  boundary  layer,  calculated  from  the  formula - *  Wc~~*  ^iel 

follows  from  (4l)  and  (44).  For  comparison,  Fig.  5  gives  the 
corresponding  curves  obtained  in  [5]  using  the  velocity  profile  that 
was  determined  semiempirically  by  Harris  [6], 


u05max 


•,  which 


. 

•  1J.7 


^8*  Limiting  values  of  surface 
:.riction  and  boundary  layer  thickness 

df  *  0,  u  *  const,  B  «  const], 

U6  r 

1  -  i  cr-  max  •  i  _  w  min  „  .  .  ,  .. 

18  v  -  ,  2  Solid  line: 

pu 

data  from  this  article;  dashed  line: 
data  from  [5]. 


FTD-HT-23-3-69 


The  presently  examined  case  corresponds,  approximately,  to  flow 
on  the  dielectric  wail  of  the  inlet  section  of  an  MHD  channel  with 
parallel  walls.  With  a  sufficiently  large  value  of  parameter  N,  the 
boundary  layers  on  the  opposite  walls  of  the  channel  do  not  join. 

In.  this  case  the  friction  across  the  channel  in  the  section  of 
developed  flow  will  be  determined,  according  to  (37),  by  the 
relationship 


,t0 

Sr 


( 


-2.  Jk-Y 

■  *  ’ 


where  h  is  half  the  height  of  the  channel;  Re^  = 

Ur\6  p 

here  — - -  and  introducing  the  designations 

we  get 


Substituting 


2  2 

=  aB  h 
pv  1 


Prom  this  formula  it  follows  that  there  is  no  friction  in  the  channel 
2  2 
when  ^  >  (Harris  [6]  obtained  >  0.6). 


For  a  more  precise  solution  of  the  problem  of  the  boundary  layer 
in  the  channel  we  must  taken  into  account  the  increased  velocity  on 
the  channel  axis.  The  velocity  along  the  axis  in  the  inlet  section 
of  the  channel 

W) 


which  follows  from  the  law  of  conservation  of  matter.  We  will  consider 
that  magnetic  induction  is  sufficiently  high  and  that  the  boundary 
layers  on  opposite  walls  of  the  channel  do  not  join.  Then  having 
substituted  into  (39)  the  velocity  from  (47),  we  get  an  expression 
for  the  friction  coefficient: 


(48) 


In  the  first  approximation  we  can  substitute 

v* 

-  which  corresponds  to  a  plate,  i.e.,  use 


into  (28)  the  value 
the  expression 


FTD-HT-23-3-69 


15 


Which  follows  from  (4l)  and  (45:).  for  the  section  of  developed  flow 
(6*  »  const)  we  have,  from  (48)  and  (4i),, 


x.„-8iv 


s£ 
'  ¥ 


~  4.Wt*™ 

“  0,525*^"  ’ 

l~  k*” 


(49) 


where  kc  ■  is  determined  from  (44)  (when  b  «  0),.  In  Pig.  6  the 

IUaX 

friction  coefficient  calculated'  from  (49)  is  compared  with  Likodis * 
experimental  results  [7]. 


.  ‘  x 

ajno 


t 

Fig.  6,  Comparison  of  theoretical 
results  for  friction  coefficient  in 
an  MHD  channel  with  experimental  data  [7]. 

1  -  Re  «  2.5*10^;  2  -  Re  =  4.5-1011. 

4  .  4 

Squares:  Re  «  2.5*10  ;  crosses:  Re  »  4.5*10  . 

The  segment  of  ppssible  comparison  is  bounded,  along  the 
abscissa,  from  below  by  the  value  of  N  for  which  <$max  «  h  and  which 
can  be  determined  approximately  ^from  Fig.  5,  and  from  above  by  the 
value  of  N  corresponding  to  the  transition  of  turbulent  flow  into 
laminar.  From  the  comparison  it  follows  that  (49)  corresponds  to 
experimental  data  within  5%. 

Let  us  also  note,  as  follows  from  (42),  that  if  the  velocity  on 
the  outer  boundary  of  the  boundary  layer  u  and  induction  B  are 
connected  by  the  relationship 
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(»+*>afe  — 


(50> 


the  solution  of  (42>  can  be  written  in  the  form 


K  h. 


To  satisfy  condition  (5,6)  we  set,  e.g,. 


Then  we  get 


>c  ' 


0J9  ReJ*+‘ 


1 


This  method  of  solving  the  problem  of  the  turbulent  boundary  layer 
with  electrical  and  magnetic  fields  is  also  applicable  during  variable 
electrical  conductivity  and  another  dependence  of  u  on  x.  In-  these 
latter  cases  the  obtained  differential  equations  should  be  solved 
numerically. 
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